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Abstract 



Local percolation probabilities are used to characterize the connectivity in 
porous and heterogeneous media. Together with local porosity distributions 
they allow to predict transport properties [HJ. While local porosity distribu- 
tions are readily obtained, measurements of the local percolation probabili- 
ties are more difficult and have not been attempted previously. First mea- 
surements of three dimensional local porosity distributions and percolation 
probabilities from a pore space reconstruction for a natural sandstone show 
that theoretical expectations and experimental results are consistent. 



According to the mixing law of local porosity theory [fl| p~0|| the effective frequency de 



pendent dielectric function e e (u) of a heterogeneous mixture may be calculated by solving 
the integral equation 
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In ([!]) the local porosity distribution n(4>) and the local percolation probability \(<p) are 
geometrical input functions which will be defined shortly. The functions £i(ou), and 62(00) 
are the frequency dependent dielectric functions of the pure materials. If ei(u),i — 1,2 and 
A(0) are known then (D) gives a prediction for the effective dielectric function s e (uj). 

Measurements of the local porosity distribution fi((p) are readily obtained from thin 



sections |6||||ri|] . Measurements of the local percolation probabilities A(0), on the other 
hand, are more difficult, and have not been reported. The main purpose in this paper is to 
present preliminary measurements of local percolation probabilities for a natural sandstone 
specimen. 

Define P to be the pore space of a two component porous medium. The complement 
of P in R 3 is called the matrix space M. More generally P and M may represent the 
two components in a heterogeneous medium. The local porosity 0(K) measured within a 
measuring region K is defined by 

where V(K) is the volume of a set K. Let g» G Z 3 denote the lattice vectors of a simple 
cubic lattice, and let 



K (r, L) = {r' G M 3 : |r- - n\ < L/2, i = 1, 2, 3} 



(5) 



denote a cube of sidelength L centered at r. 

Given the notation introduced above the local porosity distribution may be defined as 



i M 

M0, L) = lim — £ K<t> - 0(K(L gl , L))) (6) 



if the limit exists. Here S(x) is the Dirac 5-function. In practice the sample is finite, and 
hence the limiting process terminates after a finite number M of measurement cells. The 
resulting histogram is used as an approximation for //. 

Local porosity distributions quantify the fluctuations in volume fractions. To characterize 
the transport properties,however, it is necessary to quantify the degree of connectedness of 
the porous medium. Two points inside the pore space P are called connected if there exists a 
path entirely within the pore space that connects the two points. A cubic measurement cell 
Kj in a cubic partitioning is called percolating in the x-direction if there exists a path within 
P R Kj connecting those two opposite faces of Kj that are perpendicular to the x-direction. 
Percolation in the y- or ^-direction is defined analogously. The local percolation probability 
\ x {4>, L) in the x-direction is defined as the fraction of measurement cells that are percolating 
in the x-direction and have a local porosity <fi. The local percolation probabilities X y (<p,L) 
and \ z ((f>,L) are defined analogously. 

To measure /) and A(0, L) in practice it is necessary to reconstruct the three dimen- 
sional pore space P for digital processing. This was done by the method of serial sectioning 



for a Savonnier oolithic sandstone specimen [p~2],[3~3|] . This type of sandstone exhibits oomoldic 



porosity with ellipsoidal pores and grains of sizes between 200/im and 300/mi [0. After 
filling the stone with a coloured epoxy resin the specimen was cut and polished. The pol- 
ished surface was photographed before removing another layer of material parallel to first 



surface. The second surface was then polished and photographed before iterating the pro- 
cedure. A total of 99 sections was prepared. The distance between consecutive planes was 
about 10/im. Each photograph represented an area of roughly 2cmx2cm. The photographs 
were scanned and thresholded into a binary image with 1904 x 1904 pixels. This corresponds 
to a resolution within each image of roughly 10/im. Figure [I] shows a typical binary image 
from within one of the 100 planes measuring roughly 1cm along each side. 

The family of three dimensional local porosity distributions L) for this sample is 
obtained straightforwardly. The results are displayed in Figure |2] for various L. Ref. |J has 
investigated the question to what extent the three dimensional local porosity distribution 
can be measured from two dimensional sections. While this is not generally possible, it was 
found that a two dimensional measurement combined with an appropriate rescaling of the 
length L gives good results. 

A preliminary measurement of the local percolation probabilities has also been carried out 
for the oolithic sandstone specimen. The results are shown in Figure |3| for cubic measurement 
cells of size L = 32. The curve X x (4>; 32) is shown in black, the curve X y (4>; 32) is shown in 
medium gray, and the curve A z (0; 32) is shown in light gray. The local porosity distribution 
for L = 32 is reproduced as the inset. 

While the data show much scatter these preliminary estimates for the local percolation 
probabilities indicate that the sandstone has anisotropic connectivity. The same trend is 
observed for smaller as well as for larger values of L, although the strength of the anisotropy 
appears to change slightly. 

The data are in agreement with the theoretical expectation that the local percolation 
probabilities should increase monotonically from to 1. The shape of the curves conforms 
to the one expected for systems with a percolation threshold. An example is the grain 
consolidation model with random packings For small values of L, however, the shape 
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of the local percolation probabilities resembles that for the central pore model All of 
these conclusions are very preliminary. More work is necessary to corroborate them, and to 
validate the measurement. 

ACKNOWLEDGEMENT: The authors are grateful to Chr. Ostertag-Henning, E. 
Haslund, Dr. U. Mann, Prof.Dr. B. N0st, Prof.Dr. D.H. Welte and Prof.Dr. R. Koch 
for discussions, and to Norges Forskningsrad and the Forschungszentrum Jiilich for partial 
financial support. 



4 




FIG. 1. Digitized thin section image of Savonnier oolithic sandstone The pore space P is 
coloured black, the matrix space is rendered white. 
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porosity 

FIG. 3. Local percolation probability functions \ X (4>,L) (black), A y ((f),L) (medium gray), 
and X z ((f>,L) (light gray) for the Savonnier oolithic sandstone shown in Figure [l] using a simple 
cubic lattice of measurement cells with side lengths (lattice constant) L = 32. The local porosity 
distribution n(<p, L) for L = 32 is shown as the inset. 
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